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ABSTRACT 



We present detailed analyses of oxygen K absorption in the interstellar 
medium (ISM) using four high-resolution Chandra spectra towards the X-ray 
low-mass binary XTE J1817-330. The 11-25 A broadband is described with a 
simple absorption model that takes into account the pileup effect and results in 
an estimate of the hydrogen column density. The oxygen K-edge region (21-25 A) 
is fitted with the physical warmabs model, which is based on a photoionization 
model grid generated with the xstar code with the most up-to-date atomic 
database. This approach allows a benchmark of the atomic data which involves 
wavelength shifts of both the K lines and photoionization cross sections in order 
to fit the observed spectra accurately. As a result we obtain: a column density 
of JV H = 1.38 ± 0.01 x 10 21 cm" 2 ; ionization parameter of log£ = -2.70 ± 0.023; 
oxygen abundance of Aq> = 0.689to]oioi anc ^ ionization fractions of O I/O = 
0.911, O II/O = 0.077, and O III/O = 0.012 that are in good agreement with 
previous studies. Si nce the oxygen abundanc e in warmabs is given relative to 



the solar standard of iGrevesse fc Sauvall (119981 ). a rescaling with the revision by 



Asplund et al. 



(]2Q09l ) yields A Q = 0.952+q;q 2 3, a value close to solar that rein- 
forces the new standard. We identify several atomic absorption lines — Ka, K/3, 
and K7 in O I and O II; and Ka in O III, O VI, and O VII — last two probably 
residing in the neighborhood of the source rather than in the ISM. This is the first 
firm detection of oxygen K resonances with principal quantum numbers n > 2 
associated to ISM cold absorption. 



Subject headings: atomic data — atomic processes — ISM: general — X-rays: general 
— X-rays: binaries — stars: individual (XTE J1817-330) 
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Introduction 



X-ray studies of interstellar absorption have the potential to provide information that 
is not accessible by other techniques. This includes the relative abundances of a wide range 
of ion stages of interstellar cosmic elements (e.g. oxygen) and the abundances relative 
to H and He of elements with atomic number Z > 6. X-ray absorption can also provide 
signatures of the binding of these elements in molecules or solids; the inner-shell electronic 
transitions are key diagnostics since the ionization state or chemical binding shifts the line 
energies by a predictable amount. The ongoing effort to explore and exploit this effect in 
interstellar studies involves accumulation of synthetic spectra for inner-shell absorption of 
ions and neutrals, based on both atomic calculation and laboratory experiment, in addition 
to testing them against available observed interstellar X-ray absorption spectra. 

As part of the effort to accumulate the required atomic data, t he atomic database of the 



xstar photoionization modeling code ( iBautista &: Kallman 



200ll ) has been systematically 



improved in the past ten years in order to study the K lines and edges in the high-quality 
X-ray astronomical spectra obtained from the Chandra and XMM-Newton satellite-born 
observatories. This database includes K-vacancy levels, wavelengths, A-values, radiative 
and Auger widths and high -energy photoioniza tion cross sections for complete isonuclear 



sequences with Z < 30 (see 



Palmeri et al. 



2012 



and references therein). 



With this new database it has been possible, for instance, to estimate the model 
efficienc y of iron K-line emiss ion and absorption in terms of the ionization parameter and 



density (jKallman et al 



20041 ). It has been shown that the centroid of the Ka unresolved 
transition array, the K/3 energy, and the ratio of the Kai to K«2 components are useful 
diagnostics of the ionization parameter. It was also found that the many strongly damped 
resonances below the K-ionization thresholds lead to edge smearing, which can certainly 
hamper the astrophysical interpretation of the absorption features. A synthetic spectrum 
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has been modeled by lKallman et al.l (120091 ) for the black-hole X-ray transient GRO J1655-40 
which constrains the odd-Z element abundances (11 < Z < 27) and outflow parameters 
of the associated warm absorber. By analyzing the oxygen K-absorption structure in the 
interstell ar medium (I S M) fr om the XMM-Newton spectrum of the low-mass X-ray binary 



Sco X-l, 



Garcia et al. 



(120111 ) reproduced both the K edge and O I Ka absorption line, 
thus evaluating the impact of the atomic data on the interpretation of observations. In the 
present report we continue this work with a more stringent benchmark of the xstar atomic 
database, made possible by Chandra spectra of a low-mass binary displaying imprints of 
the ISM oxygen K absorption with higher resolution. 



Juett et al. 



( 120041 ) measured high-resolution spectra of the interstellar oxygen K-shell 
absorption edge from seven X-ray binaries — Cygnus X-l, Cygnus X-2, 4U 1636-53, 
4U 1735-44, GX 9+9, 4U 1543-624, and 4U 1820-30— using the Chandra High Energy 
Transmission Grating Spectrometer (HETGS). The oxygen column density Nq was therein 



calculated using the optical depth at 21.7 A and 
assuming the oxygen abundance from I Wilms et al. 



; he cro ss section at this wavelength 



(120001 ) . With this value, the hydrogen 



column density iVn was then determined. The ls-2p transitions of O II and O III were 
identified, finding ISM abundances for O II and O III relative to O I of ~ 0.1 and < 0.1, 
respectively. This became the first estimate of the O II/O I and O III/O I abundance ratios 
in the ISM. 



By observing the oxygen K-shell edge toward galaxy clusters, 



Baumgartner fc Mushotzky 



( 120061 ) showed that the ISM hydrogen X-ray column density was in close agreement with 
the 21 cm radio value for columns less than approximately 0.5 x 10 21 cm" 2 . For higher 
column densities, the former was higher by as much as a factor of 2.5, indicating substantial 
absorption besides that due to neutral hydrogen (probably from clouds of molecular 
hydrogen). An average ISM oxygen abundance of 0.99 solar was found relative to the 
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photospheric value (8.66 ± 0.05) by lAsplund et al.l (120041 ) which implied a high gas to dust 
ratio. 



Pinto et al.l (120101 ) reported an spectral model suggesting that, along the line of sight 
towards an X-ray source (GS 1826-238), the ISM was more complex than just simple neutral 
gas. Their gas model consisted of three components: cold gas (T « 5.8—10 x 10 3 K); warm 
ionized gas (T « 1-6 x 10 4 K); and hot ionized gas (T « 2 x 10 6 K). It was also shown 
that the column densities Nr of these three components span two orders of magnitude, 
where the cold ga s contributed ~ 90—95% , the warm gas ~ 5—10%, and the hot gas ~ 1% 



to the total A^ ot 



Costantini et al 



(120 12l ) analyzed the ISM using the bright X-ray binary 
4U 1820-30 source, concentrating on the mildly ionized cold gas where the model consisted 
of mildly ionized gas, dust, and mildly ionized outflowing gas. In the oxygen edge region, 
dust compounds were included which did not show features below 23.7 A. Finally, oxygen 
was found to be overabundant by a factor of 1.23 solar. 

In order to model the oxygen K absorption features in the ISM with xstar, we need a 



Chand ra spectrum with high counts. In the catalog of low-mass X-ray binaries by 



Liu et al. 



(120071 ). we found in XTE J1817-330 a suitable source with an average c ount number of 
« 8 x 10 6 counts/s. Discovered by the Rossi X-ray Timing Experiment (jRemillard et al. 



20061 ). XTE J 18 17-330 is a low- mass X-ray binary displaying a primary component which 



has an X-ray sp ectrum anc 



as a black hole (ISala et al.l 



varia bility behavior which are consistent with its identification 
20071 ). Its location, I = 359.81 and b = -7.99 in Galactic 



coordinates, implies a line of si ght in the Galactic ce nter direction, a region that dominates 



the total mass of the Galaxy (jStasiriska et al. 



20121 ) providing a good opportunity to 
study the ISM absorption f eatures. This sourc e has been obse rved in different spect ral 



ranges , namely the 



radio (IRupen et al. 



optic al ( jTorres et 



al.l 120061). near in 



20061 ). and X-ray (IRov et al. 



201 lh . 



rared 



D'Avanzo et al. 



Sala et al. 



2006|) 



(120071 ) fitted its X-ray 
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spectrum finding an absorbing hydrogen column of 1.55 x 10 21 cm 2 and identifying the 
following oxygen absorption lines: O I Ka( 23.52 A); O II Ka/23.3 5 A); O III Ka(23.13 A); 



Gierlihski et al. 



O I K/3(22.91 A); and O VII K/3(21.60 A), 
column density of iVfi = 1-08 ± 0.05 x 10 21 cm -2 , and 



Roy et al. 



( 20091) estim ated a hydrogen 



(120111 ) investigated the 
timing characteristics in the spectra using a column of iVfi = 1.2 x 10 21 cm -2 . 

The outline of the present report is as follows. In Section 2, we describe the xstar 
warmabs physical model which is used to obtain the final spectrum fits, in particular the 
details concerning its atomic database. The data reduction aspects of the observations, 
e.g. the pileup effect, are reviewed in Section 3 while the results of the broadband fits and 
resolved absorption features are given in Sections 4-5. Section 6 is dedicated to comparisons 
of the hydrogen column density that allow an estimate of its reliability. Section 7 describes 
the oxygen column densities derived from our fits. Finally, a discussion and conclusions are 
presented in Section 8. 



2. XSTAR/ warmabs 

xstaiCj is a code package designed to determine the physical conditions in partially 
ionized gases. It calculates the temperature, ionization, excitation, absorption, and emission 
due to neutral and all ionized species of elements with atomic number Z < 30. As many 
relevant physical processes as possible are taken into account, assuming a time stationary 
balance among them and a Maxwellian electron velocity distribution. Important processes 
for many problems relevant to X-ray astronomy include photoionization, electron impact 
collisional ionization and excitation, and radiative and dielectronic re combination. More 



complete descriptions of xstar and its atomic database are given by 



Kallman fc Bautista 



1 http:/ /heasarc. nasa.gov/lheasoft/xstar/xstar. html 
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feOQlh and 



Bautista &: Kallmanl (120011 ) . In this work we employ the warmabs utility 



which allows for the calculation of xstar models within the widely used X-ray spectral 
fitting package xspecJ^. It utilizes precalculated tabulations of the ion fractions and level 
populations in order to calculate synthetic emission and transmission spectra. It assumes 
that the gas responsible for emission or absorption has a uniform ionization and temperature 
throughout, although it is possible to superimpose two or more components and simulate a 
non-uniform situation. 



2.1. Atomic data 



The K photoabsorption data sets for oxygen ions with electron numbers 2 < N < 8 



Garcia et al. 



( 120051 ). Energies 



contained in the xstar atomic database were generated by 
for both valence and K-vacancy fine-structure levels, transition wavelengths, ^/-values, 
and radiative and Auger widths were computed with the atomic stru cture package s 



autostructure ( Eissner et al. 



1974 ; 



Badnell 



1986 



19971 ) and hfr (ICowanl Il98lh . 



Electron correlation effects were taken into account by means of configuration-interaction 
expansions within the n < 3 complexes, and relativistic corrections were included with a 
Breit-Pauli Hamiltonian. The accuracy of the theoretical level energies and wavelengths 
does not generally match that attained in experiments (see Table [T|), and the scanty 
availability of measurements for the K shell of O ions with N > 3 limits the implementation 
of term-energy corrections and wavelength shifts that can be introduced in order to 



compensate for the theoretical shortcomings. Furthe rmore, as also s 



experimental misassignments, such as that quoted by 



Gu et al. 



rown in Table [IJ 



(120051 ) for O VI, are not 



uncommon. (The discrepancy in O III is believed to be a typo.) 



2 http:/ /heasarc. gsfc.nasa.gov/xanadu/xspec/ 
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High-energy photoionization cros s sections in the K-edge region for oxygen ions 
were o btained by 



Garcia et al. 



codes ( IScott fc Burke 



1980 



( 200 5T) with bprm, the Breit-Pau 



Scott &: Taylor 



1982 



Berrington et al. 



i i?-matrix suite o 



1987 



Scat on 



19871 ). A 



conspicuous effect in the resonance structure of these cross sections that must be treated 
in detail is Auger damping. It is the result of the dominant spectator KLL Auger decay 
channels of a photoexcited K-vacancy state, say, 



[ls]2p At np -> 2p^ 2 np + e" 

->• [2s]2p^~ 1 np + e" 
->• [2s] 2 2p M np + e" 



(1) 
(2) 
(3) 



over the participator Auger KLn channels 



[ls]2p^np -> 2p^ 1 +e" (4) 
->• [2s]2p" + e- , (5) 

which leads to edge smearing by resonances with symmetric profiles of nearly constant 



width. They are handled with i n the -R-matrix package 



( iGorczyca fc Robicheaux 



1999: 



Gorczyca fc McLaughlin 



jy me ans of an optical potential 



20001 ) . Again, the accuracy of the 



resonance energy positions is limited by the absence of measurements which would enable 
empirical adjustments of the many K-vacancy target thresholds (i.e. the series limits) in 
the close-coupling expansion. 



3. Observations and data reduction 



In Tabled we list the specifications of the four observations of XTE J1817-330 obtained 
by Chandra using the High Energy Transmission Grating Spectrometer (HETGS) in 
combination with the Advanced CCD Imaging Spectrometer (ACIS). The oxygen edge is 
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accessible with the HETGS exclusively via the Medium Energy Gratings (MEG). Given 
the large flux of this source (F = 1350 /iJ, integrated in the 2-10 keV energy range), 
these observations were taken in continuous clocking mode which significantly increases the 



20081 ). The pileup 



temporal resolution in order to minimize the pileup effect (jCackett et al.l 
is an inherent feature of CCD detectors, such as those of the ACIS instrument, which occurs 
when two or more photons are detecte d as a single event possibly causing a deformation in 



the level and shape of the continuum (iMiller et al. 



2006b). It is stronger in the MEG than 



in the HEG due to the lower dispersion and higher effective area of the formecl. Although 
the pileup is usually not an issue near the oxygen edge, it may be for the continuum which 
is used to establish the hydrogen column density by fitting the 11-25 A interval; hence, we 
consider the pileup effect in this range by applying the correction model simple_gpile2 . si 



(jHanke et al. 



20091 ) before spectral fitting. Among the most popular computer packages 
for fitting X-ray data are xspec and ISIS, and we have used the latter (version 1.6. 43) to 
include the pileup model. 

Even though the spatial resolution is reduced to one dimension when the continuous 
clocking mode is used, there is essentially no difference with the timed exposure mode in 
the background extraction or in the data reduction process. Hence, we have reduced the 
data sets using the standard CIAO thread^. In some cases, the zero-order data were not 
telemetered; for these observations the zero-order position was then estimated by finding 
the intersection of the grating arms using the f indzo algorithrd_\ 

Figure [1] shows the light curves of the four observations ObsID 6615, 6616, 6617, and 



3 http: / / cxc.harvard.edu/ ciao/ahelp/acis_pileup.html 
4 http:/ /space. mit.edu/cxc/isis/ 
5 http: / / cxc. harvard. edu/ciao4.4/threads/gspec. html 
6 http: / / space.mit.edu/ cxc/ analysis/findzo/ 
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6618 in units of counts/s. It may be seen that the average counts/s is approximately 
constant for all, and therefore, we use the averaged spectra for this work. We have carried 
out the following spectral fitting procedure: firstly, each spectrum was rebinned with a 
0.1 A bin in order to fit the continuum in the 11-25 A wavelength range so as to find the 
hydrogen column density using the TBnew model; then, the default spectral resolution was 
restored, and the oxygen edge region (21-24 A) was analyzed with the warmabs model 
freezing the column density at the value obtained in the preceding step. 



Broadband fit 



Following 



Yao fc Wangi ( 20061 ). we rebinned each spectrum with a 0.1 A bin size 



throughout the spectral range to decrease the number of absorption lines and improve 
the broadband fit. We fitted the four observations simultaneously in the interval 11-25 A 
using the simplegpile2 (TBnew (powerlaw) ) model, where simplegpile2 is a convolution 
model to account for pileup effects and TBnew is an X-ray absorption model that includes 



abundances for elements from H to Ni ( Wilms et al. 



20001 ) . Also we added Ga ussians to fit 



the re maining absorption line s after rebinning. W e used abundances specified in 



Wilms et al. 



( 120001 ) and cross sections by 



Verner et al. 



(119961 ). The O, Ne, and Fe abundances were 



handled as free parameters, and for the analysis \ 2 statistics was employed. We fixed 
the absorption parameters in all the observations, and varied the power-law, pileup, and 
Gaussian parameters for each case. 

Figure |2]shows the fit in the interval 11-25 A using the simplegpile2 (TBnew (powerlaw) ) 
model. Although the fit is carried out simultaneously for the four observations, they 
are plotted separately to improve clarity. Panels (a), (b), (c), and (d) correspond to 
observations ObsID 6615, 6616, 6617, and 6618, respectively (see Table [2]). In each panel, 
black data points correspond to the observation while the solid red lines represent the 
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best-fit model for each case, and the base plots show the fit residuals in units of x 2 - in eac h 
spectrum the K edges of Ne and O are respectively located at ~ 14.3 A and ~ 23.2 A. The 
O I and O II ls-2p absorption lines at 23.5 A and 23.35 A are also clearly observed. We 
improved the fit by including Gaussians to model the outstanding absorption lines: Ne IX 
Ka (» 13.43 A); Ne III Ka (» 14.50 A); Ne II Ka (» 14.60 A); O VIII Lyct (» 18.95 A); 
O VII Ka (« 21.58 A); and O II Ka (» 23.34 A). Residuals indicate that the model fits 
the data satisfactorily with deviations near the Ne edge (~ 14.3 A) due to background (i.e. 
clocking the whole chip) or an unaccounted hot pixel. We also find high residuals in the 
oxygen absorption region; in particular, the O I Ka and O II Ka are not well modeled. 

Fit results for the individual free abundances, column density, power law, and pileup 
are presented in Table [3j The fit statistics give a reduced chi-square of \ 2 — 1-208. 
Differences in the power- law parameters are due to the coun t numbers in each observation. 
We obtain abundances for O, Ne, and Fe higher than solar (jWilms et al.ll2000i ). especially 
for Ne (factor of 2), and a column density of Ah = 1.66ig^ x 10 21 cm_2 which is somewhat 
larger than those estimated in other spectral fits and by measurements of the 21 cm line. 
This will be further discussed in Section [61 

The nonlinear pileup convolution model simple_gpile2 . si exponentially reduces the 
predicted count rate Rj(X) (in units of counts/s/A) according to 

R(\) = R(\)exp(-PR bat (\)) (6) 

where R to t{X) is the total spectral count rate and j3 is treated as a fit parameter. The model 
is only applicable to first order grating spectra although it employs information from the 
second and third order grating spectra[_. With this model we obtained an average (3 = 0.050 
according to which the spectra have a pileup degree greater than 25% for 14 < A < 16 A 



7 http: / / cxc.harvard.edu/ciao/ahelp/acis_pileup.html 
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with the highest value of p = 42% at fh 14.3 A for the observation with the largest number 
of counts (ObsID 6615). This value shows the need to include the pileup effect in the 
11-25 A region to obtain a good fit. However, for all the observations, we have a pileup 
degree lower than 5% in the 21-24 A wavelength range which corresponds to the oxygen 
absorption region and, therefore, its effects were ignored. 



5. Absorption features 



In order to identify all the absorption features present in the 21-25 A wavelength 
region, we first applied a functional model consisting of a power law for the continuum and 
Gaussian profiles to describe the absorption lines. Figure |3] shows the spectral fit of the four 
Chandra MEG observations of XTE J1817-330 in the oxygen absorption region. Although 
shown in separate panels for clarity, the fit was performed simultaneously for the four 
observations. Nevertheless, the power-law parameters (photon index and normalization) 
are allowed to vary independently for each data set. Several absorption features are clearly 
observed, the most outstanding being the O I and O II Ka lines. The base plot in each 
panel shows in units of x 2 the fit residuals as a histogram. High residuals are found around 
24.8 A, and the model underestimates the number of photons around 23.2 A. 

Table H] shows a list of the lines obtained with the functional fit including the centroid 
wavelengths derived from the Gaussi a n pro file parameters. It also lists the wavelengths 



for the lines measured by 



SalaetaL 



( 120071 ). Note that our wavelengths are within the 
error bars of the values reported by these authors, which are much larger given the lower 
spectral resolution of the XMM-Newton instruments. We obtain a reduced chi-square 



of x 2 



1.434. The O 



declared of ISM origin (jSala et al. 



Ka is found at 23.502 ± 0.001 A which has been previously 



20071 ). Furthermore, we have also detected four other 



O I absorption lines corresponding to resonances with higher principal quantum number n: 
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two pairs of K/3 (n = 3) and K7 (n = 4) resonances associated to the 4 P and 2 P K-hole 
core states. These are found at 22.884 ± 0.004 A, 22.790 ± 0.001 A, 22.686 ± 0.004 A, and 
22.609 ± 0.007 A, respectively. In this set, only the first K/3 (belonging to the 4 P core state) 
ha s been prev i ously detected at 22.91 ± 0.03 A in the XMM-Newton observation analyzed 



by 



SalaetaL 



( 120071 ). This is then the first firm detection of oxygen K resonances with 



n > 2 associated to ISM cold absorption. 



The line at 23.358 ± 0.002 A is also of ISM origin flJuett et all 120041 ) and corresponds 
to O II Ka. As in the case of O I, resonances with n > 2 were also detected for O II. 
These were found at 22.280 ± 0.003 A and 22.101 ± 0.005 A corresponding to K/3 and K7, 



respectively. Although the O III Ka resonance is a triplet, only two absorption 



Sala et al. 



eatures 



(120071 ) only 



were detected at 23.104 ± 0.005 A and 23.054 ± 0.001 A while 

found one line at 23.13 ± 0.09 A identified as O III Ka. Finally, we have also located Ka 
transitions for both O VI (22.022 ± 0.003 A) and O VII (21.589 ± .003 A), from which only 



the strongest of the two (O VII) was reported by 



Sala et al 



(|20o3) at 21.609 ± 0.06 A. 



The functional fit described above illustrates the complexity of oxygen K absorption. 
There are (at least) two different regimes: a cold component of mostly neutral gas and a hot 
component with highly ionized gas. Previous studies have suggested a multiple-phase ISM; 
however, it is not clear if the two high-ionization lines from O VI and O VII are intrinsic 
to the source or in fact arise in the ISM. Thus we direct our study to the cold phase. 
Furthermore, since cold absorption is not exclusively due to neutral oxygen, this means 
that a model such as TBnew will not completely represent the observed features. We have 
then simultaneously fitted the four exposures in the 21-25 A wavelength range using the 
powerlaw*warmabs physical model. The advantage of using warmabs is that it contains the 
most recent atomic data for all the ions in the oxygen isonuclear sequence, thus enabling 
the fit of all the absorption lines from O I, O II, and O III. While implementing this model, 
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the column density is held fixed at the value obtained in the broadband simultaneous fit 
(iVn = 1.66 x 10 21 cm -2 ). The ionization parameter and the oxygen abundance are treated 
as free parameters but the same for all observations, while the abundances of all the other 
chemical elements are held fixed at their solar values. The quantities describing the power 
law (photon index and normalization) are also taken as free parameters but independently 
for each observation. 

Figure H] shows the results of this fit. As before, we present for clarity the spectra 
with the best-fit model in separate panels although the fit is performed simultaneously. 
Since most of the relevant absorption features in the spectra are due to neutral, singly, and 
doubly ionized O ions, the ionization parameter that describes the fit tends to be somewhat 
low: log £ = —2.729. In this case, the warmabs model does not show any absorption due to 
O VII; therefore, this line is still represented with a Gaussian profile. We have found large 
residuals around most of the oxygen absorption features, the most prominent being around 
the O I and O II Ka lines. The fit statistics gives a poor reduced chi-square of \ 2 — 1-769 
which forced us to revise the atomic database. 

A close examination of the resonance positions from the photoabsorption cross sections 
used in warmabs in comparison with those obtained from the spectral fits with Gaussian 
profiles reveals the inaccuracies of the atomic data. Figure [5^ shows with solid colored 



lines the atomic cross sections used in warmabs for O I, O II, and O III (IGarcia et al 



20051 ). Vertical dashed lines are placed at the positions of the absorption features found 
with the functional fit described above (see Table H]). The comparison shows that not 
only the positions of the O I and O II Ka resonances are displaced with respect to the 
observed data but also those with n > 2. The s olid black line repr esents the experimental 



O I photoabsorption cross section measured by IStolte et al 



(119971 1 (unfortunately, there 



are no laboratory measurements available for the other species) which displays a similar 
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wavelength shift with respect to the observed reson ance positions. This wavelength offset 
has been previously discussed by 



Juett et al. 



( 120041 ) based on the Ka position observed in 



the Chandra spectra from several sources, all of which match the position we have found 
using Gaussian profiles. For O I, the comparison of the theoretical cross section with the 
laboratory measurement shows that, except the Ka peak, all the other resonances agree 
very well. This implies that the relative position of the Ka with respect to the other 
resonances also needs to be adjusted. This is also the case for O II since the wavelength 
difference of the theoretical and observed Ka positions is not the same as those for the 
rest of the resonances of the same ion. Consequently, we have decided on the one hand to 
adjust the positions of the Ka resonances (which are treated as lines in warmabs), and on 
the other, to shift the whole cross sections for both species so as to obtain the best possible 
agreement with the observed lines. The new positions for the O I and O II Ka lines are 
23.502 A and 23.343 A while the cross sections were shifted towards shorter wavelengths 
by 0.033 A and 0.079 A, respectively. No correction was applied to the O III curve since 
the observed lines are weaker, in addition to the fact that in this case the Ka resonance 
is a triplet. Figure |5b shows the atomic absorption cross sections after the wavelength 
corrections are applied. In other words, the Chan dra observed li ne po sitions suggest not 



Garcia et al. 



(120051 ) must be adjusted 



only that the theoretical resonance positions from 
but a lso that the wavelength scale of the only available experiment for O I (iStolte et al 
19971 ) needs to be shifted by 33 mA. 



Possible effects due to the instrument wavelength calibration, model uncertainties, 
and Doppler shifts due to motion of the gas in the ISM are important for accurate 
line identification. The instrumental resolution for the first-order spectra from MEG is 
AA ~ 23 mA, which is better than the largest shift r equire d for our theoretical cross 



sections (33 mA). In fact, according to 



Garcia et al. 



((20051), the uncertainties in the 



theoretical cross sections for O ions with electron occupancies N < 4 were estimated to be 
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~ 50 mA, com p arab 
of 



Stolte et al 



(IJuett et al 



e with the wavelength uncertainties in the laboratory measurements 
1997f ). Assuming a velocity dispersion for the ISM of < 200 km s _1 
20041 ). the largest wavelength shift due to Doppler effects for the O I Ka is 
AA = 14 mA, which is smaller than the instrumental resolution. Hence, we can safely 
ignore the effects of ISM gas motion in the light of sight. 

Figure [6] shows the fit results with warmabs after the theoretical line positions have been 
corrected, where the residuals are now significantly reduced and more evenly distributed. 
This is reflected in better fit statistics (reduced x 2 — 1-245) and a well-constrained ionization 
parameter (log£ = —2.699 ± 0.023) which is consistent with the low-ionization oxygen lines 
observed. As before, the O VII Ka resonance at 21.589 A is not well represented with this 
model and must be fitted with a Gaussian profile. We have also noticed that the absorption 
feature at 22.022 A is not entirely due to the O III K5 resonance and is thus fitted with a 
Gaussian profile, confirming a blend of both O III K5 and O VI Ka. In this final analysis, 
we also treated the column density value as free parameter in order to increase the accuracy 
of the fit. We find a column density of Nn = 1.38 ± 0.1 x 10 21 cm -2 that is lower than 
the value previously mentioned in connection with the broadband fit (see Section 4). The 
oxygen abundance obtained with the warmabs model is Ao = 0.689^g'o^o relative to solar, 
also lower than the value for the best fit with TBnew (Ao = 1.178 ± 0.22). All the best-fit 
parameters of the warmabs are summarized in Table [5j 



6. Hydrogen column density 

Table [6] shows a comparison of the hydrogen column density obtained in this work 
with previous estimates including measurements using the 21 cm line. By means of the 
simultaneous broadband fit described in Section HI we have obtained a column density of 
N H = 1.661 

aol x cm 2 wn i cn is somewhat larger than those quoted in previous spectral 
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fits. iMiller et all ()2006aJ) reported a hydrogen column density of iVn = 8.8—9.7 x 10 20 cm 2 
using a 50 ksec Chandra spectrum, a value smaller than those obtained from Galactic 



1 I surveys, e.g. 1.5 8 x 10 cm (jDickey fc Lockman 



( iKalberla et al. 



optical spectrum, 



1990|) and 1.39 x 10 21 cm 



20051). Using t h e equ ivalent widths of several interstellar bands in the 



Torres et al 



( 120061 ) fo und the column de nsity 



iV H = 1-3 x 10 21 cm" 2 . Moreover, both iRvkoff et al. 



;o be in the ran ge 



(120071 ) and 



Gierlihski et al. 



(I2008h 



gathered acceptable fits of several Swift XR 



of N R = 1.2 x 10 21 cm" 2 . 



Gierlihski et al 



spe ctra with a fixed Galactic column density 



(120091 ) used the same XRT spectra combined 



with data from UVOT (The Ultraviolet and Optical Telescope) to derive a column density 
of ]Vh = 1.08 ± 0.05 x 10 21 cm -2 from broadband spectral fits. 



The agreement between TBnew and 



Sctlct 6t ell. 



(120071 ) is around 7% despite the fact that 



their fits include data from XMM-Newton EPIC-Pn (0.6-10.0 keV), RGS1 (0.3-2.0 keV), 
and OM covering a much wider energy range; fur thermore, they mo del the foreground 



absorption with TBabs, an older version of TBnew (jWilms et al 



20001 ). On the other hand, 



the present fits of the oxygen K region (21—25 A) with the warmabs physical model yield a 
lower column density (Nh = 1-38 ±0.1 x 10 21 cm" 2 ) than that derived from our broadband 
fit and TBnew. This discrepancy is expected due to the different atomic data sets involved, 
particularly since TBnew only includes the photoabsorption cross section for neutral oxygen. 
Taking into account the scatter of the previous column densities in the literature and the 
fact that we have adequately modeled the absorption features from O I, O II, and O III, we 
would expect the value derived from warmabs to be more reliable than that from TBnew. 



7. Oxygen column density 



In order to derive the O column density along the line of sight of XTE J18 17-330, 
we have calculated the curve of growth for each of the oxygen Ka transitions using the 
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atomic data from iGarcfa et al.l (120051 ) . and equivalent widths (EWs) have been obtained 
using the functional fit described in Section Table [7] shows a co mparison be t ween the 



present Ka EWs for O I, O II, O HI, O V 



XTE J1817-330. EWs listed in 
as well as those by 



Yao et al. 



Juett et al- 



and Q VII and those by 



Sala et al 



( 120071 ) for 



(120041 ) for seven other sources are also included 



is around the lower limits of ISala et al. 



(2009Q for Cygn us X- 2. For O I Kg, ou r EW (51 ± 5 mA) 



Juett et al. 



(200 



Juett et al. 



), an d for O II Ka, 



(120041 ) associated with 



(120071 ) and 

the present value of 52 ± 5 mA is within the range that 
the ISM. However, we find the Ka lines in both O I and O II to be saturated which could 
influence our derived EWs. Furthermo re, the Ka EWs in O III, O VI, and O VII a re in good 



agreem ent with the values obtained by 



Juett et al. 



(120091 ). 



(120041 ) . 



Sala et al. 



(120071 ) . and 



Yao et al. 



Table El shows a comparison of the oxygen column densities obtained from the EWs and 
the warmabs model fit. As previously mentioned in Section [51 for the ionization parameter 
found in our best warmabs fit, only O I, O II, and O III are included in the model. As 



a refe rence, we have also included the oxygen column density for O VI from I Savage et al. 



( 120031 ). obtained from FUSE data for four sources located near the same Galactic latitude as 
XTE J1817-330. The curve-of-growth values for O I, O II, and O III result in a total column 
density of No = Aoi + Aon + Nom = 5.71 x 10 17 cm -2 , in relatively good agreement with 
that derived from the ionic fractions in the warmabs fit, namely Nq = 6.41 x 10 17 cm -2 . 
However, the ratios O I/O II are significantly different: 1.4 and 11.9, respectively. Due to 
line saturation, the measured EWs are likely to be underestimated, in particular for the 
O I Ka line. This leads to a higher uncertainty in the column and ion fractions estimated 
from the EW measurements. On the other hand, since the parameters in the warmabs fit 
are constrained by all the lines as well as the K edge, we are more confident about their 
reliability. In the case of the O VI and O VII Ka transitions, because the EWs are on 
the flat section of the curve of growth, they strongly depends on velocity dispersion. For 
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O VI Ka, the oxygen column density spans from No vl = 71.12 ± 29.15 x 10 15 cm -2 to 
^Ovi = 5-06 ± 2.53 x 10 15 cm -2 when using dispersion velocities of v = 20 km s _1 and 
v = 200 km s _1 , respectively. The smalles t valu e is still about one order of magnitude 



larger than that reported by 



Savage et al. 



(120031 ). making difficult to constrain the the 



velocity dispersion. Also, no real improvement is achieved by considering larger dispersions, 
as the column density is reduced very little once the velocities are larger than 200 km s -1 . 
This may be an indication that the observed O VI line is likely to originate from the 
neighborhood of XTE J1817-330 rather than from the ISM. In the case of O VII, the 
oxygen column density varies from Nq vii = 4.48 ± 1.84 x 10 18 cm -2 (v = 20 km s^ 1 ) to 



Nr 



4.5 ± 1.8 x 10 ib cm" 2 v = 200 km s 



8. Discussion and conclusions 



There are several interesting conclusions emerging from this work that configure a set 
of useful guidelines in the analysis of ISM absorption features in Chandra and XMM-Newton 
spectra. In order to resolve such spectral signatures, bright sources with high count rates 
are desirable. In the present XTE J1817-330 study case, we were fortunate to find four 
good quality spectra that were fitted simultaneously so as to obtain improved statistics. 
However, it was realized that the pileup effect must be taken into account in order to obtain 
reasonably reliable hydrogen column densities. This effect was treated adequately with the 
nonlinear pileup convolution model simple_gpile2 . si, whereby the pileup was estimated 
by fitting the continuum to be on average around 25%. The resulting hydrogen column 
density is JV H = x 10 21 cm 2 , whi ch is somewhat (< 7%) higher than those by 



Kalberla et al. 



(120051 ) and 



SalaetaL 



(120071 ). the former obtained from 21 cm data. 



The absorption features in the oxygen region (21-25 A) were first picked up with a 
functional model (power law and Gaussian profiles) which yielded Ka lines of O I, O II, 



-21 - 



O III, O VI, and O VII previously observed by ISala et al.l (120071 ). but for the first time, 
also K resonances of O I and O II with principal quantum number n > 2. These line 
identifications at least confirm a two-phase plasma with both neutral (mostly) and highly 
ionized components. In the present work we have been mainly concerned with the former. 
In this respect, since neutral oxygen is not the sole charge state responsible for the cold 
oxygen absorption, the four exposures are then simultaneously fitted with a representative 
physical model referred to as war mabs, which c ontai ns the atomic data for the oxygen 



isonuclear sequence computed by 



Garcia et al 



(120051 ) . In this step, the hydrogen column 
density is held fixed at the previously determined value of Nn = 1.66 x 10 21 cm" 2 while 
the ionization parameters and oxygen abundance are treated as free parameters. A detailed 
comparison of this model with observations brings forth the inaccuracies of the atomic 
data, which led us to introduce wavelengths shifts in both the Ka line positions and 
photoabsorption cross sections of O I and O II. 

The scope of this corrective procedure is severely limited by the general poor 
availability of laboratory data for the oxygen inner-shell, but it certainly leads to 
improved fits and statistics; in particular, to a well-constrained ionization parameter of 
log£ = —2.699 ± 0.023 consistent with the absorption lines observed. In order to further 
improve the fit, the column density is allowed to vary resulting in a more reliable value 
(Nfi = 1.38 ± 0.1 x 1 21 cm -2 ) which is within 13% of previous estimates (IKalberla et al. 



2005 



Salaetal 



20071 ). Furthermore, the fit yields an oxygen abundance of Aq = 0.689 



+0.015 
-0.010 



and ionization fractions of O I/O = 0.911, O II/O = 0.077 and O III/O = 0.012. Regarding 
the oxygen abundance, it may be clarified that the quoted v alue is given relative to 



the solar standard (8.83 ± .06) of 



Asplund et al 



Grevesse fc Sauvall (119981 ). If it is rescaled to the 



(120091) of 8.69 ± 0.05, then an ISM abundance 



recently revised standard by 
(Aq = .952ln m°) close to solar is obtaine d that is in excellent accord with the averaged 



value by 



Baumgartner &: Mushotzkyl (120061 ). and can be included as further support for the 
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solar abundanc e revision. 
those found by lJuett et al 



Moreov er, our ionization fractions are in good agreement with 
fl2004T ) of O II/O I « 0.1 and O III/O I < 0.1. 



The functional fit allows us to derive oxygen column densities by means of a curve of 
growth method using the EWs for each of the observed Ka absorption lines, which can 
then be compared with the values obtained from the warmabs fit. Even though we find by 
both methods a relatively consistent total column density (5.71 x 10 17 cm -2 in the curve of 
growth and N = 6.41 x 10 17 cm -2 in the warmabs model fit), the O I/O II ratio for the Ka 
lines is significantly lower (an order of magnitude) in the former. This is attributed to the 
high saturation level in both lines which leads to underestimates of the measured EWs, and 
we are therefore confident that the results from the warmabs model fit are more reliable. 
In the case of O VI, the width of the line strongly depends on the velocity dispersion of 
the gas. For the velocity dispersion values expected in the ISM (y ~ 20 — 200 km s _1 ), 
the EWs obtain ed with our function al fit yield oxygen column densities far above from the 



FUSE values by 



Savage et al 



(120031 ). This may be an indication that the O VI observed 
absorption in the XTE J1817-330 spectra is bound to occur near the source rather than in 
the ISM. 

We find acceptable fits to the observed oxygen K-edge spectrum with a model consisting 
of atomic oxygen ions if the energies are slightly adjusted. However, it is also of great 
interest to detect or set limits on features due to other phases of interstellar oxygen and 
other elements, namely molecules or solids (dust). It is expected that these phases comprise 
a significant fraction of the total oxygen column in the ISM. Evidence for this includes 
observed depletions of atomic oxygen along m any lines of sight relative to certain elements 



( iJenkins 



2009 



Yao et al. 



2009 



Whittet 



20101 ). together with observations of molecular 



features in other wavelength bands. An example of a potentially detectable signature of such 
material is the K resonance of the CO molecule expected at a rest wavelength of 23.21 A 
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( IBarrus et al.lll979l ). This line is separated from the nearest feature we clearly detect by 
approximately 100 mA, which can be compared with the magnitude of the wavelength 
shifts we adopt in order to get acceptable fits, the largest of which is 33 mA. Therefore, we 
conclude that the existing atomic absorption data are not sufficiently uncertain to mask the 
presence of this molecular signature. The search for molecular and solid features in X-ray 
spectra remains of great interest, both as tests of the atomic data and of the models for the 
ISM. 

We believe that the present benchmark has produced an improved version of the 
warmabs model which can now be used for a more extensive and reliable study of oxygen 
K photoabsorption towards other X-ray sources. This will be the next step of the current 
study where the oxygen abundance variations in the ISM are certainly of interest. In a 
similar fashion, we also intend to evaluate the Ne and Fe edge regions for ongoing physical 
model refinement. 
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2012 to the European Space Astronomy Centre (ESAC), Madrid, Spain, the latter funded 
by the COSPAR Fellowship Program. Warm hospitality and tutelage at these institutions 
are kindly acknowledged, in particular from Andy Pollock and Carlos Gabriel at ESAC. We 
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Fig. 1. — Light curves of the four observations, ObsID 6615, 6616, 6617, and 6618, in units 
of counts/s. The approximately constant average counts/s that prevails is an indication of 
a low variability degree. 
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Fig. 2. — Chandra MEG spectra of the X-ray binary XTE J1817-330 simultaneously fitted 
in the 11-25 A region using the TBnew model, (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 
6617. (d) ObsID 6618. See Table [2]for further details pertaining these observations. In each 
panel, the black data points are the observations while the solid red lines correspond to the 
best-fit models. Residuals are shown as histograms for each case in units of x 2 - Note the 
large residuals in the oxygen K-edge region, particularly in the neighborhood of the O I and 
O II Ka lines at « 23.50 A and 23.35 A, respectively. 
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Fig. 3. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a simple functional model (power law and several Gaussian 
profiles), (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. The observed 
absorption features are interpreted as: the neutral O K-edge (~ 23.1 A); the Ka, K/3, and 
K7 lines of O I (w 23.50 A, w 22.88 A, and w 22.68 A, respectively); the Ka, K/3, and K7 
lines of O II (w 23.35 A, w 22.28 A, and « 22.10 A, respectively); and the Ka lines of O III, 
O VI, and O VII (» 23.10 A, w 22.02 A, and w 21.58 A, respectively). See Table H for more 
details. 
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Fig. 4. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a powerlaw*warmabs physical model, (a) ObsID 6615. 
(b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 
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Fig. 5. — Comparison of the theoretica 
(green), and O III (blue) computed by 



photoabsorp t ion c ross sections for O I (red), O II 



Garcia et al 



(j2005l ) which are i mplemented in the 



Stoltc et al. 



(1997) 



warmabs model. The black solid line is the laboratory measurement by 
showing peaks of O I and O2. The vertical dashed lines are placed at the wavelengths of the 
observed absorption lines (see Table H]). Panel (a) displays the original cross sections while 
panel (b) shows the same curves after the wavelength shifts are applied to both O I and O II 
(see text for details). 
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Fig. 6. — Spectral fit of the Chandra MEG observations of XTE J1817-330 in the oxygen 
absorption region (21-25 A) using a powerlaw*warmabs physical model, after the wavelength 
corrections in the O I and O II atomic cross sections are introduced, (a) ObsID 6615. (b) 
ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 
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Table 1. Comparison of EBIT and theoretical wavelengths (A) for K lines 


T 

ion 


T 1 1 / J\ 

Lower level [J) 


Upper level (J ) 


BB11 


1 heory 


VI 


ls 2 2sd/2) 


ls2s2nfl/2 3/2) 


22.0194 ± 0.0016° 


22.00 c 








22 374 ± 008 b 


22.05 d 










22.03 6 


V 


ls 2 2s 2 f0) 


1s2s 2 2d(1) 


22 374 ± 003 a 


22.33 c 








22.370 ±0.010 6 


22.35 d 










22.37 e 


IV 


ls 2 2s 2 2p(l/2,3/2) 


ls2s 2 2p 2 (l/2,3/2) 


22.741 ± 0.005 6 


22.78 c 










22.73 d 










22.75 e 


III 


ls 2 2s 2 2p 2 (l,2) 


ls2s 2 2p 3 (l) 


22.071 ± 0.006 6 


23.08 c 










23.05 d 










23.07 6 



i EBIT measurement (ISchmidt et al.l 120041 ) 



b EBIT measurement (jGu et al 



2005|) 



hullac calculation (IBehar fc Netzerll2002l ) 



d i?-matrix calculation (IPradhan et al 



20031 ) 



2 hfr calculation ( iGarcfa et al. 



20051 ) 
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Table 2. Chandra observations used in this paper 



ObsID 


Date 


Exposure (ks) 


Instrument 


Grating 


Label 


6615 


2006 Feb 13 


18 


HETG-ACIS 


MEG 


Observation 1 


6616 


2006 Feb 24 


29 


HETG-ACIS 


MEG 


Observation 2 


6617 


2006 Mar 15 


47 


HETG-ACIS 


MEG 


Observation 3 


6618 


2006 May 22 


51 


HETG-ACIS 


MEG 


Observation 4 
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Table 3. Broadband simultaneous fit parameters 



Model Parameter Value 



TBNew 


Nh (10 21 cm' 2 ) 


1 


fifi +0.03 
DD -0.04 




TBNew 


Neon abundance (Anj e ) 


2 


258 


±0 


080 


TBNew 


Iron abundance (Ap e ) 


1 


302 


±0 


050 


TBNew 


Oxygen abundance (Aq) 


1 


178 


±0 


022 


Power-law 


Normalization 


9 


416 


±0 


064 a 






7 


039 


±0 


041 6 






4 


613 


±0 


042 c 






1 


481 


±0 


001 d 


Power-law 


Photon Index 


1 


545 


±0 


030° 






1 


618 


±0 


030 6 






1 


735 


±0 


029 c 






2 


242 


±0 


038 d 


SimpleGpilc2 


P 





050 


±0 


001" 









051 


±0 


001 b 









050 


±0 


001 c 









049 


±0 


001 d 


Reduced chi-square 


x 2 


1 


208 







Note. — Abundances relative to the solar values of 



Wilms et al 



(2000) 



a ObsID 6615 
b ObsID 6616 
c ObsID 6617 
d ObsID 6618 
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Table 4. Absorption line assignments 



Ion 






Ttp n <>it inn 


A a (A) 




VII 


Is 2 


1 S- 


- [ls]2p 4 P° 


21.589 ±0.003 


21.60 ±0.06 


VI 


2s 2 


S - 


[ls]2s2p 2 P° 


22.022 ±0.003 




on 


2p 3 


4 S° 


- [ls]2p 3 4p 4 P 


22.101 ±0.005 




on 


2p 3 


4 S° 


- [ls]2p 3 3p 4 P 


22.280 ± 0.003 




01 


2p 4 


3p 


- [ls]2p 4 ( 2 P)4p 3 P° 


22.609 ±0.007 




01 


2p 4 


3p 


- [ls]2p 4 ( 2 P)3p 3 P° 


22.686 ±0.004 




01 


2p 4 


3p 


- [ls]2p 4 ( 4 P)4p 3 P° 


22.790 ±0.001 




01 


2p 4 


3p 


- [ls]2p 4 ( 4 P)3p 3 P° 


22.884 ±0.004 


22.91 ±0.03 


III 


2p 2 


3p 


- [ls]2p 3 3 P° 


23.054 ±0.001 


23.13 ±0.09 


III 


2p 2 


3p 


- [ls]2p 3 3 D°, 3 S° 


23.104 ±0.005 


23.13 ±0.09 


on 


2p 3 


4 S° 


- [ls]2p 4 4 P 


23.358 ±0.002 


23.35 ±0.03 


01 


2p 4 


3p 


- [ls]2p 5 3 P° 


23.502 ±0.001 


23.52 ±0.02 



Present work 



ISflla et all fl2007h 
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Table 5. warmabs simultaneous fit parameters 



Model 


Parameter 




Value 


Warmabs 


7V H (10 21 cm" 2 ) 


1 


38 ±0.01 


Warmabs 


Log ionization parameter (log £) 




2.699 ±0.023 


Warmabs 


Oxygen abundance (Aq) 





fi o Q +0.015 
uoa -0.010 


Power-law 


Normalization 


2 


384° 






2 


544^ 






1 


511 c 









685 d 


Power-law 


Photon Index 


3 


044° 






2 


479 6 






2 


720 c 






2 


617 d 


Reduced chi-square 


x 2 


1 


245 



Note. 



O xygen abundance relative to the solar value of 



Grevesse fc Sauvall a998l) 
a ObsID 6615 
b ObsID 6616 
c ObsID 6617 
d ObsID 6618 
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Table 6. Hydrogen column density comparison 



Method 



N H (10 21 cm" 2 ) 



TBnew fit a 


1 66+ 03 
± ' DD -0.04 


warmabs fit" 


1.38 ±0.1 


21 cm survey 6 


1.58 


21 cm survey c 


1.39 


Spectral model fit^ 


1.55 ±0.05 


Spectral model fit e 


1.2 


Spectral model fit-f 


0.88 - 0.97 


Spectral model fit 9 


1.0 - 3.0 


Spectral model fit' 1 


1.08 ±0.05 



^Present work 



TOickev fc Lockmar] ( 1990) 



Kalberla et al 



(2005) 



d Sala et al 



1 20071 1 



c 


Rvkoff et al. 




2008 


). 


Rov et al. 



(2007) 



Gierlinski et al 



[ 201ll l 



f Miller et al 



i 2006al l 



? Torres et al 



(2006) 



h Gicrliiiski ct al 



(2009) 
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Table 7. Equivalent widths comparison for the oxygen Ka transitions in the ISM 



Source 


O I Ka 


O II Ka 


O III Ka 


O VI Ka 


O VII Ka 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


XTE J1817-330 a 


51 ±5 


52 ±5 


16 ±5 


14 ±5 


54 ±4 


XTE J1817-330 6 


70 ±20 


40 ±20 


30 ±30 




40 ±30 


Cyg X-2 c 


67±11 


26±17 


9±8 






Cyg X-2 d 


50 ±3 


28 ±6 




11 ±2 




4U 1543-624 c 


68±11 


41±11 


18±5 






4U 1820-30 


70±20 


41±17 


20±10 






4U 1735-44 c 


80±30 


40±20 


24±10 






GX 9+9 c 


80±30 


50±19 


27±13 






4U 1636-53 c 


90±30 


60±20 


33±13 






Cyg X-l (ObsID 3407) c 


95±17 


20±11 


35±15 







a Present work 



11 Sala etal 



(2007) 



Juett et al 



(2001) 



c Yao et al 



(2009) 
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Table 8. Oxygen column density comparison 



Source 


1 


O II 


O III 


O VI 


O VII 




(10 17 cm- 2 ) 


(10 17 cm" 2 ) 


(10 17 cm- 2 ) 


(10 15 cm- 2 ) 


(10 15 cm- 2 ) 


XTE J1817-330 a 


3.18 ±0.79 


2.18 ±0.54 


0.35 ±0.08 


5.06 ± 2.53 d 


45 ± 18 d 










5.99 ± 1.9 e 


988 ± 405 e 










71.12 ± 29.15 f 


4484 ± 1838 f 


XTE J1817-330 b 


5.85 ± 1.75 


0.49 ±0.12 


0.07 ±0.02 






PG 1302-102 c 








0.16 ±0.03 




Mrk 1383 c 








0.38 ±0.07 




ESO 141-G55 c 








0.31 ±0.06 




PKS 2005-489 c 








0.60 ±0.01 





a Calculated from the EWs and the curves of growth (see section 0). 
b Derivated from the warmabs model fit using the solar value of J 



Grevesse &: Sauv al (1998) (see section 



Savage et al 



c Obtained from FUSE data by 
d Using a velocity dispersion v = 200 km/s 
c Using a velocity dispersion v = 100 km/s 
f Using a velocity dispersion v = 20 km/s 



(2003) 



